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^SJ ' There are a number of reasons to entertain the possibihty that locaHty is violated 

on microscopic scales, for example through the presence of an infinite series 
of higher derivatives in the fundamental equations of motion. This type of 
nonlocality leads to improved UV behaviour, novel cosmological dynamics and is 
O ■ a generic prediction of string theory. On the other hand, fundamentally nonlocal 

ZJ [ models are fraught with complications, including instabilities and complications 

rG ^ i in setting up the initial value problem. We study the structure of the initial value 

problem in an interesting class of nonlocal models. We advocate a novel new 

CNJ ■ formulation wherein the Cauchy surface is "smeared out" over the underlying 

<^ I scale of nonlocality, so that the the usual notion of initial data at t = is 

^^ • replaced with an "initial function" defined over — A/^^ < t < where M is 

the underlying scale of nonlocality. Focusing on some specific examples from 

rs^ ■ string theory and cosmology, we show that this mathematical re-formulation 

has surprising implications for the well-known stability problem. For D-brane 

^\ • decay in a linear dilaton background, we are able to show that the unstable 

(^ [ directions in phase space cannot be accessed starting from a physically sensible 

initial function. Previous examples of unstable solutions in this model therefore 
correspond to unphysical initial conditions, an observation which is obfuscated 
in the old formulation of the initial value problem. We also discuss implication 
of this approach for nonlocal cosmological models. 
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1. Introduction 

Nearly all dynamical systems encountered in physics are local. Time evo- 
lution is usually described by differential equations involving at most two time 
derivatives, and the t > dynamics depend only on a small amount of data 
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specifying the state of the system at i = 0. There are, however, a number of 
reasons to entertain the possibihty that locahty is violated on microscopic scales, 
for example through the appearance of an infinite series of higher derivatives 
in the fundamental equations of motion. This type of nonlocality is a generic 
feature of many particle physics models including string theory [1| , field theories 
living on noncommutative space-time (2| and also models with a minimal length 
scale [SJ. In fact, it has been argued that any Quantum Field Theory (QFT) of 
gravity must be nonlocal at the most fundamental level [J]. 

The underlying nonlocality of string theory manifests itself strikingly through 
the infinite derivative structure of String Field Theory Q (SFT) and also toy 
models such as the p-adic string [y, |7| or strings quantized on a random lattice 
[8[. These constructions make it possible to study both nonperturbative and 
off-shell string physics, and have therefore played a key role in describing the 
decay of unstable D-brane configurations (tachyon condensation) . The dynam- 
ics of tachyon condensation [9| are relevant for both formal string theory and 
also cosmology, because brane/anti-brane annihilation [lOJ marks the endpoint 
of popular D-brane inflation models jll| . 

Aside from their "bottom up" motivations, field theories with infinitel y rn an^j 
derivatives are interesting due to improved Ultra- Violet (UV) behaviour lli 
and novel cosmological dynamics [li, [H, [11, [13, E, [H, S, [2l|, [H, [H . In [2l[ 



it was shown that nonlocal theories can support inflation even on a very steep 
potential, alleviating the fine tuning problems that plague inflationary model 
building. This mechanism leads to a robust prediction for large nongaussian- 
ity in the Cosmic Microwave Background (CMB) which will be falsifiable (or 
verifiable!) with Planck [H,!!!!- 

In spite of these motivations, physical applications of nonlocal field theories 
are fraught with complications. Perhaps the most serious concerns are the 
following. 

1. Stability: If the equation of motion admits more than two initial data 
then the "extra" degrees of freedom can be interpreted as physical ex- 
citations which carry wrong-sign kinetic energy and render the Hamilto- 
nian unbounded from below [J, |l2|, |2J, |25|. The quantum theory con- 
tains ghosts, while the classical theory is plagued by Ostrogradski (higher 
derivative) instabilities. 

2. Predictivity: If the equation of motion requires infinitely many initial 
data, then there is an independent concern. By suitable choice of the 
infinite free parameters of the solution it might be possible to arrange 
for nearly any time dependence over an arbitrarily long interval. In this 
case the initial value problem would be completely bereft of predictivity. 
(Throughout this paper we understand the word "predictivity" in this 
sense.) 

Clearly, very much is at stake. There is evidence that both of these problems 
might manifest themselves in string theory [J], for example through the wild 
oscillatory dynamics of tachyon condensation in flat space [7| (however, see 
0). 



As emphasized in [2J] , the presence of ghost- hke instabihties in a nonlocal 
QFT is intimately related to the structure of the initial value problem for the 
classical equations of motion. In this work, we consider a new formulation of 
the initial value problem for a broad class of nonlocal theories which incorporate 



dissipation [27|, |28|, |29| into the dynamics (as is unavoidable in a cosmological 



context). We show that, in such models, the infinite tower of initial conditions 
d^ (p a.t t = can be swapped for the freedom to fix the time evolution over 
some interval At = T by requiring 

cP{t)^Mt) on te[-T,0] (1) 

with 00 (i) the "initial function" and T~^ related to the energy scale at which 
nonlocal effects become important. This new formulation has considerable in- 
tuitive appeal since it makes manifest the nonlocal structure of the theory; the 
Cauchy surface has been "smeared out" over a finite interval set by the under- 
lying scale of nonlocality. 

This mathematical re-formulation of the initial value problem suggests a 
rather novel new perspective on the key problems of predictivity and stability 
discussed above. Although the initial function is arbitrary from a purely mathe- 
matical viewpoint, we point out that (/'o(i) rnust be subjected to some auxiliary 
constraints in order to ensure that the — T < t < evolution is physically sen- 
sible. For example, in the case of D-brane decay (poit) should be restricted to 
the physically meaningful region of the potential. The physical interpretation of 
these constraints is completely obfuscated in the "old" formulation of the initial 
value problem wherein one fixes infinitely many derivatives at t = 0. 

Remarkably, these mild physical restrictions on 0o (t) have the effect of con- 
straining the theory to a stable subset of its solution space. Although the 
Hamiltonian is globally unbounded below, the unstable directions in phase space 
cannot be accessed from an allowed initial configuration. This restriction is non- 
perturbative and retains information about the full nonlocal structure of the 
equations. Moreover, the constrained theory makes robust predictions for the 
macroscopic {At ^ T) dynamics. 

In light of these results, we argue that there exists a broad class of stable, 
predictive nonlocal models whose equations of motion admit infinitely many 
initial data. A specific example is provided by D-brane decay in the background 
of a linear dilaton [271, |28|, |29| . Previous examples of unstable solutions in this 
model were implicitly assuming unphysical initial conditions. This perspective 
on the nature of the instability was not appreciated in previous studies because 
the initial value problem was formulated in a highly unnatural way. 

Our results help to establish the theoretical consistency of nonlocal cosmo- 
logical models, and also provides a plausible mechanism by which the consistent 
inclusion of closed strings might enable SET to evade the Ostrogradski instabil- 
ity. (This has been conjectured as a possible solution of the instability problems 
that plague SET, however, to our knowledge no concrete mechanism has previ- 
ously been proposed.) Since we do not rely on any intrinsically stringy properties 
of the theory, it will be straightforward to apply our approach to more general 



nonlocal models, in cosmology or elsewhere. At the practical level, our analysis 
provides a new tool-box of efficient and numerically stable techniques for solving 
nonlinear equations with infinitely many derivatives. 

The organization of this paper is as follows. In section [5] we discuss some 
general properties of delayed differential equations which will be important for 
the subsequent analysis. In section |3] we re-visit models of light-like D-brane 
decay in light of the formulation ([ij of the initial value problem, showing how 
macroscopic predictivity and stability can be salvaged. In section|3]we show how 
the simple-minded inclusion of dissipation can lead to delay-type nonlocality in 
a very general nonlocal QFT framework and we show how the formulation ^ of 
the initial value problem can be implemented. In section[5]we study a particular 
model, cosmological D-brane decay, showing again how stability/predictivity 
can be salvaged. Finally, in section [6l we conclude. 

2. Nonlocal Models with Delays 

In this section we will provide a brief discussion of some general features 
of models whose nonlocality arises purely through delay operators of the form 
g-TOt ^g .^g .^jjj ggg shortly, this type of nonlocality is often associated with 
the presence of some source of dissipation in the dynamics (such as cosmological 
friction). 

2.1. The Method of Steps 

Consider a prototype model whose dynamics are governed by the following 
delay differential equation 

N 

J2a,4''U{t)^Fm^T)] (2) 

n=0 

Here F [x] is some nonlinear forcing term and we assume that T > 0. Equation 
(m is a nonlinear differential equation of infinite order, which is easy to see by 

noting the identity (j){t - T) = e-'^^^(t>{t) = Y.n -^^[^^("VlO- Infinite order 
equations have been studied in some detail using the formal generatrix calculus 



23, and also I25|,|3l|. (See [32, 33, 34, ^] for different approaches.) 



The model ([2]) belongs to a special class of nonlocal theories, referred to 
as having "compact support" in |4], which can be formulated in a surprisingly 
simple and intuitive manner using the method of steps (see also !29J). The key 
observation is that equation ^ defines a mapping from functions on the interval 
[t — T, t] to the interval [i, t + T]. We can therefore break the solution (/)(i) into 
segments 4>i{t) defined over contiguous intervals of length T. Explicitly, we set 

(f,{t) ^ (j),{t) on te[{i-l)T,iT] (3) 

where the integer i = 0, 1, 2, • • • labels the segment under consideration. This 
procedure for splitting a solution up into segments is illustrated schematically 



in Fig. [TJ The dynamical equation (I5|) provides a recursion relation to obtain 
the i-th solution segment: 



N 



Y, andrUt) = F [</.,_i(t -T)], z = 1, 2, . 



(4) 



n=0 



On the first segment {i — 0) we have an arbitrary "initial function" <j)Q{t). More 
on this shortly. 




Figure 1: Schematic illustration of the procedure for splitting a function up into solu- 
tion segments defined over contiguous intervals of length At = T. 



Starting from a given 4>o{t), the stepping procedure (JH) is very easy to im- 
plement. Each solution segment is obtained by solving an inhomogeneous linear 
differential equation whose source term is constructed from the previous solu- 
tion segment. At the points t — iT where the segments are joined we should 
impose boundary conditions that (j)(t) be continuous and as smooth as possible. 
Since the differential operator on the left-hand-side of dH) is A^-th order, we can 
demand that 



9t^"V.+i(t) 



- ;^(") 



d^Ut) 



t=iT 



for n = 0, • • • , iV - 1 



(5) 



This matching should be performed at each knot: i = 0, 1,--- The method 
of steps is nonperturbative and solutions obtained using the recursion (JH) are 
often much more efficient and numerically stable than those obtained using 
alternative approaches (such as higher order perturbation theory [241] or the 
diffusion equation technique |18|). 

Notice that, in general, the 7V-th derivative of (j){t) will be discontinuous at 
t = 0. This discontinuity will propagate into higher derivatives in subsequent 
segments. The only way to avoid this effect is by choosing (poit) to be a solution 
of ([2]). In practice, however, we can minimize derivative discontinuities by using 
perturbation theory to construct the initial function. None of our qualitative 
results in this work will rely on the smoothness properties of the solutions in 
any way. 



2.2. The Initial Value Problem 

The initial value problem associated with 1^ can be formulated as 

(j){t) = Mt) on te[-T,0] (6) 

where the initial function (j)o{t) specifies the state of the system on the first 
segment. Equation 1^ can be regarded as "smearing out" the Cauchy surface 
over a finite width. Notice that ^ is not the unique formulation of the initial 
value problem for this model. If we have adopted perturbation theory (instead of 
the method of steps) it would have been most natural to specify the derivatives 
d\ (j) at i = 0, following [24|. Although both views of the initial value problem 
are mathematically valid, the formulation ([6]) is much more natural here because 
it makes the manifest the delay-type nonlocality of the underlying theory. 

The initial value problem (jH]) is unusual. The conventional notion of causal- 
itjo is violated, because the evolution depends on past history. Suppose that 
we are observing the dynamics of (t>{t) over a time scale short compared to the 
delay. At ^ T. Equations ([2]) and dH) imply that the system will be directly 
infiuenced by events which occurred very far in the past, before the start of the 
"experiment" . No finite amount of information about the state of the system 
at t = is sufficient to uniquely determine the i > evolution. Indeed, equa- 
tion ([5]) means that we could arrange for any time dependence over an interval 
At = T, regardless of the specific details of the dynamical equation ^. Clearly, 
the nonlocality of the system is associated with a breakdown of predictivity over 
microscopic scales, At <C T. 

On the other hand, if we observe the system for a macroscopic time scale. 
At ^ T, then the situation is different. We retain a weakened version of causal- 
ity because the "memory" is small as compared to the dynamical scales of 
interest. We could hope to measure the evolution over an entire segment and 
use this as an initial function to robustly predict the subsequent dynamics. (In 
principle this still requires making infinitely many measurements and, in a real- 
istic model, it might only be possible to reconstruct the initial function to some 
limited accuracy. Any uncertainty in 0o(i) propagates into all future predic- 
tions. In practice, this effect is usually tiny because the late-time dynamics are 
insensitive to the details of the initial function.) 

Notice that, although the initial function may be ambiguous in some cir- 
cumstances, it is not completely arbitrary. We should subject (jioit) to some 
auxiliary constraints in order to ensure that the —T < t < evolution was 
sensible. The specific nature of these constraints will, obviously, depend on the 
system that we are modeling. Notice that we have already seen an example of 
this necessity: arbitrary choices of 4>o{t) will lead to derivative discontinuities. 
We will see shortly how the imposition of physical restrictions on 0o (t) can have 
very significant implications for the t > dynamics of the system. 



^For our purposes, the conventional notion of causality can be caricatured as "the present 
determines the future" . 



3. Delayed Dynamics in String Theory 

In the last section we introduced a general class of models with delay-type 
nonlocality, equation ([2]) • We argued that ^ provides the most natural formu- 
lation of the initial value problem. Let us now consider two explicit examples 
from string theory and show how this new perspective on the initial value prob- 
lem provides invaluable insight into the key problems - stability and predictivity 
- discussed in the introduction. 

3.1. An Example from p-adic String Theory 

Let us now consider an explicit example of the type ^ whose nonlinear 
dynamics are particularly transparent. Our example is based on p-adic string 
theory [6|, a toy model of the bosonic string which is known to reproduce a 
number of key features of the full string theory. (See [36| for a review of some 



of these features and [3a, |37[ for a discussion of the p-adic theory at finite 
temperature and implications for the cosmological constant problem.) The p- 
adic theory contains a single scalar field (pix), the open string tachyon, that 
encodes the instability of a space-filling D-brane. Following [27| we adopt light- 
cone coordinates x^ = -75 (2;" ± x^) so that the metric takes the form ds'^ = 
— (da:")^ -|- (dx^)"^ = —2dx^dx~' and assume a linear dilaton profile: ^{x) = 
—V^x~. In this closed string background the p-adic Lagrangian is J27| 



L = 



9 



■W^^^^r' 



p) 



where p is a prime number, gp is related to the open string coupling, we have 
set a' = m^'^ = 1. The potential is 

gle-^*^-V{<t>) = +\^^-^^r^' (8) 

which is plotted in Fig. [2j If we assume a light-like tachyon profile </> — (t>{x'^) 
then the equation of motion derived from ^ takes a remarkably simple form 

0(x+)=0P(x+-T) (9) 

where the delay is 

r = y+lnp (10) 

Equation © admits two constant solutions: (/) = and cj) = 1 (also (j> = —1 
when p is odd) corresponding, respectively, to the local maxima and minima 
of the potential. Physically, the solution (/> = 1 corresponds to state with an 
unstable space-filling D-brane while (f) = is the empty state containing neither 
D-brane nor open string excitationso The dynamical process of rolling the 



■^This configuration is sometimes called the "vacuum" although, properly, the theory (0 
has no minimum energy state. 
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Figure 2: The potential for the p-adic tachyon field (j> taking p = 3 for illustration. 
The maxima </!) = ±1 correspond to the physical state containing an unstable D-brane 
while the minimum at = is the "vacuum" containing neither D-brane nor open 
string excitations. The unbounded regions are thought to be associated with the closed 
string tachyon. The solid part of the potential, — 1 < (/> < -1-1, denotes the range of 
field values which can be accessed by physically meaningful initial functions. 



tachyon from the unstable maximum to the minimum of the potential - tachyon 
condensation - provides a time-dependent description of D-brane decay in the 
p-adic string theory. 

The initial value problem ([6]) requires that we choose the evolution of the 
p-adic tachyon over a time interval Ax+ = T = V'^ In p. Which choices of 
(j)o{x^) will lead to a sensible physical interpretation in terms of brane decay? 
To answer this question, first notice that the potential ^ is unbounded below 
for (/) > 1 (also (j) < — 1 for p odd). This is usually thought to be associated with 
the presence of the closed string tachyon [39| and solutions rolling down the 
unbounded region are regarded as unphysical. Second, recall that near = 
there are no open string excitations (because there is no brane for them to end 
on) and hence the physical meaning of the open string variable is obscure. 
We should therefore expect that the theory ([7]) only admits a robust physical 
interpretation in terms of brane decay when the initial function satisfies: 



< Mx^) < 1 



(11) 



When p is odd we can continue this to —1 < (j)o{x^) < 0. Thus, only the solid 
region of the potential curve plotted in Fig. [5] can be accessed by physically 
meaningful initial functions. 

Equation © can be easily solved using the method of steps. However, it 
happens that we can write down a simple exact solution which is valid on all 
segments, as was pointed out in [28[. That solution is 



<p{x ) — exp 



^/y^ 



F{x+) 



(12) 



where F{x'^) is a periodic function satisfying 

F{x+) = F{x+ - T) (13) 

(The solution (|12p is always positive, however, when p is odd we can obtain a 
second solution simply by flipping the sign.) The initial value problem ^ fixes 
the free function F{x'^) on the first segment as 

F{x+) = -e-"^/^^ In [Mx^)] on t G [-T, 0] (14) 

The periodicity relation (jl3|) then uniquely determines F{x'^) on all subsequent 
segments. 

Armed with the exact closed-form solution (J12p we can address the nonlinear 
dynamics of ([7]) is a remarkably simple way. If the initial function satisfies (|llll 
then we must have F{x'^) > and the solution (|12p approaches the minimum 
(f> = very rapidly for t >^T . This behaviour is illustrated in the left panel 
of Fig. [21 Notice the peculiar dynamics: 4'{x'^) can roll to higher potential 
energy during the rolling, even though the dynamics are governed by gradient 
flow. Such motion is possible because the kinetic energy in this theory can be 
negative. Nevertheless, the nonlinear dynamics are stable in the sense that the 
field does eventually settle down to the minimum of the potential at late times. 
Note also that these late-time asymptotics are largely insensitive to details of 
Mx+). 

The Hamiltonian associated with ([7]) is unbounded from below [28|. How- 
ever, we have just argued that the nonlinear solutions of this theory do not 
display any unstable behaviour. How can this be? To imderstand what is going 
on, we should imagine relaxing the constraint (jlip on the initial function. Now 
the dynamics associated with ([9]) are very different. If (J)q{x'^) > 1 anywhere on 
x^ e [— T, 0] then we must have F{x^) < at one or more points and conse- 
quently the solution (J12l) undergoes wild oscillations with divergent amplitude. 
This is a manifestation of the Ostrogradski instability. These dynamics are il- 
lustrated in the right panel of Fig. [31 As can be seen, the qualitative late-time 
behaviour is now highly sensitive to the details of (j)o{x^). 

In the formulation ([6]) it is quite apparent what is going wrong with the 
solutions plotted in the right panel of Fig. [3l We have allowed the tachyon to 
sample unphysical regions of the potential on the initial segment —T < t < 
0. However, this conclusion would be completely obfuscated if we adopted a 
perturbative approach. Writing 4>{x~^) = 1 + d<p{x'^) we obtain 

p^''^+ -p\S(p{x+)^0 (15) 

to linear order in S(p. The solution of (|T5|) is 

S6(x+)= V a„e^^"^^ 



Y. ane''"'-^^ (16) 
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Figure 3: The time evolution of the p-adic tachyon field for several representative 
choices of the initial function. In the left panel we consider only examples which satisfy 
the constraint < (fio{x'^) < 1. At late times the field settles down to the minimum 
of the potential and the macroscopic (Ax'^ S> T) dynamics are largely insensitive to 
the details of (j>o{x^). In the right panel we consider two examples where the initial 
function violates the physical constraint < (j}o{x'^) < 1- Now the field undergoes 
erratic oscillatory behaviour at late times and the dynamics become highly sensitive 
to the details of (j)o{x^). In both panels we set p = 3 and V"^ = 1 for illustration. The 
horizontal yellow lines illustrate the critical values tf) — and 0=1. 
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which depends on infinitely many free parameters a„. The frequencies M„ 
correspond to the zeroes of the generatix /(s) = p^ ^ — P- Exphcitly, we have 



Mn = — 



27rm 
1 



Inp 



n = 0,±l,±2,--- (17) 



The spectrum (|17p can be interpreted as a single growing mode (the usual 
tachyon) along with an infinite tower of modes which oscillate with ever-growing 
amplitudecl Half of those spurious modes jiarry wrong-sign kinetic energy and 
represent ghost-like excitations [28]. In [28] it was noted that there exists an 
"island of stability" in initial condition space. By carefully tuning the coeffi- 
cients a„ in the solution (fT6|) we can have sensible late-time dynamics. This 
tuning looks highly contrived in the perturbative approach, however, in light 
of the method of steps we see that it amounts to simply satisfying the physical 
bound (fTTj) over a microscopic time Aa;+ — T. 

Before moving on, it is worth commenting on a special class of initial func- 
tions which we have excluded from our discussion up to now. In the special 
case where (J)q{x'^) = 1 for some points on a;+ G [— T, 0], then sustained oscilla- 
tions are possible. At late times the tachyon interpolates periodically between 
the maximum and minimum, as is illustrated in Fig. |4l Taken at face value, 
this solution describes a brane which repeatedly dissolves into tachyon matter 
and then re-assembles itself. We are tempted to interpret this as an array of 
space-like brane (S-brane) |38| in the p-adic string theory. For the time being, 
however, we remain agnostic as to whether the upper limit (J)q{x^) = 1 should 
be included in the physical constraint pT|) . 

The oscillatory behaviour in Fig.|4]might seem very surprising. The dynam- 
ics of the theory ^ is purely dissipative. How can this be reconciled with the 
observation that solutions need not limit to (/> = at late times? The key to this 
puzzle is noting that the dilaton profile $(a:) is nearly constant on very small 
scales Aa;+ <^T ^ V^ . Deep in the UV the dissipative effects associated with 
the dilaton gradient must be negligible. This observation is related to the fact 
that the spectrum (|17p contains frequencies M„ which are arbitrarily large as 
compared to the energy scale associated with the damping, T^^. 

3.2. An Example from String Field Theory 

In this subsection we consider a somewhat more realistic example: the level 
zero truncation of bosonic SFT. In order to obtain a delayed formulation we 



■^The solution I I16II is completely consistent with I I12I I. This is easily seen by expanding 
Fix'^) in a Fourier series as 

F(x+) = ao+J2 «" ™s ( -—^ ) + 2^ f*.. sin ( ——- ] (18) 

When 1 2> |a,i|,|bn[ we can expand l|12|l in a Taylor series and trivially recover </'(x"'") = 
1 + E„"ne"""^. 
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Figure 4: The time evolution of the p-adic tachyon field for the special case where 
the initial function crosses (^o(a:''') = 1 at some point on a;"*" £ [— T, 0]. In this special 
case the tachyon undergoes sustained periodic oscillations between the minimum and 
maximum of the potential. It may be possible to interpret this solution as an array of 
space-like branes (S-branes). We set p = 3 and V^ — 1 for illustration. The horizontal 
yellow lines illustrate the critical values tf) — and (f} = 1. 



again assume a light-like dilaton profile of the form $(a:) = —V^x 
Lagrangian is [27] 



C = 



-lm' + l^'-^{c^°^y 



The 



(19) 



where gs is the open string coupling, the dimensionless constani|f|is C = 3\/3/4 = 
1.23 and the scalar field 4>{x), again, corresponds to the open string tachyon. 
The potential is 

(20) 



5?e-^^--l^(^) 



^^2 , ^ ,3 



which is plotted in Fig. [SJ Assuming a light-like profile cf) ~ 0(a;^) the equation 
of motion takes the form 



where the delay is 



{V+d+ - 1) (j){x+) + C^(j)^{x+ - T) = 



r = 2y+lnC 



(21) 



(22) 



Equation ([^11) cannot be solved analytically. However, the method of steps is 
easy to implement numerically. As in the p-adic case, we should subject <po{x~^) 



*Our C is related to the constant K in employed by Hellerman and Schnabl in |27| as 
C = 1/K^-^- and to the constant K employed by Beaujean and Moeller in 12911 as C = ii'B.M. 
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Figure 5: The potential for the SFT tachyon field. The maximum (j) = correspond to 
the physical state containing an unstable D-brane while the minimum at (j) — C~^ is 
the "vacuum" containing neither D-brane nor open string excitations. The unbounded 
region is thought to be associated with the closed string tachyon. The solid part of 
the potential, < </> < C~'^ , denotes the region which can be accessed by physically 
sensible initial functions. 



to some restrictions to ensure that the solutions admit a meaningful physical 
interpretation in terms of D-brane decay. What class of initial functions are 
allowed? As in our previous example, the maximum = corresponds to the 
physical configuration v^rith an unstable space- filling D-brane vi^hile the minimum 
(j) = C~^ corresponds to the empty configuration which contains neither D-brane 
nor open string excitations. At the minimum, the meaning of the open string 
variable (j) is obscure. Similarly, solutions rolling down the unbounded (/> < 
side of the potential are associated with the closed string tachyon and considered 
unphysical [39[. All told, the theory (|19p admits a robust physical interpretation 
in terms of D-brane decay only when the initial function satisfies 



< (t>o{x+) < C 



-3 



(23) 



which corresponds to the solid region of the potential curve plotted in Fig. [5j 
In addition to (P5|) . we should also ensure that the tachyon rolls down the 
"correct" side of the potential, towards the open string vacuum at (f> = C~^. 
We will discuss this point in more detail shortly, however, a sufficient (but not 
necessary) condition is the following 



0o(O) > M~T) 



(24) 



The condition (|24p allows for a variety of non-monotonic behaviours on the 
interval a;+ G [— T, 0], however, it imposes that the "net" motion of the tachyon 
over this time is to increase. As was the case with (P5)) . the condition (1^^ is not 
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onerous: one expects it to be satisfied for any solution that admits a physical 
interpretation in terms of D-brane decay. 

We have studied the nonlinear solutions obtained using the method of steps 
and found that the dynamics agree very well with the results of 28| using the 
diffusion equation approach. As long as the initial function (J3q{x~^) satisfies the 
physical constraints (|23I24[) then all solutions will eventually settle down to the 
minimum of the potential in a time comparable to the delay. This behaviour is 
illustrated in the left panel of Fig. [H] for three representative choices of initial 
function (provided explicitly in Appendix A). It should be emphasized that 
our numerical evidence for this result is much more exhausitive than the three 
specific examples plotted in Fig. [51 we have repeated this analysis for a very 
broad array of initial functions and found qualitatively similar behaviour in all 
cases. 

The Hamiltonian associated with (|19p is unbounded below, however, the 
nonlinear solutions do not display any unstable behaviour as long as (1^51) is 
satisfied. Notice also that the late-time dynamics are insensitive to the details 
of 00 (a;"*") ■ 




?-i.o 




2 4 6 



10 12 14 16 18 



Figure 6: The time evolution of the light-like tachyon field for several representative 
choices of the initial function. In the left panel we consider only examples which satisfy 
the constraint < (jf)o(x^) < C~^. At late times the field settles down to the minimum 
of the potential and the macroscopic (Ax"'" 3> T) dynamics are largely insensitive to 
the details of (j)(}{x'^). In the right panel we consider three examples where the initial 
function violates the physical constraint (l>o{x'^) > 0. The solutions may grow without 
bound and the late-time dynamics are highly sensitive to the details of 4>o{x^). In both 
panels we have set V^ — 1 for illustration and the horizontal yellow lines illustrate the 
critical values ip — and (j) ~ C~^ . The explicit choices of (l)o{x'^), for both panels, 
are provided in Appendix A. 

Just like in our last example, the dynamics are very different if we relax the 
physical constraint ([23l) on the initial function. If it happens that 0o(x"'") < 
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on a;+ S [— T, 0] then the possibility arises for (j){x'^) to roll down the unbounded 
part of the potential and grow limitlesslyO In this case the late-time dynamics 
will depend very sensitively on the details of (j>o{x^). This is illustrated in 
the right panel of Fig. [5] for three representative choices of (J)q(x^) (provided 
explicitly in Appendix A). 

Because we lack an exact analytic solution of (HJ), it is difficult to prove 
rigorously that all solutions will relax to the minimum, once (|23I24I) are satisfied. 
We have provided numerical evidence of this fact. To reinforce our conclusions, 
we now provide also analytical arguments. Firstly, we note that as soon as 
4>{x'^) approaches the basin of the closed string vacuum, (j) = C~^, then its fate 
is sealed. All solutions in the vicinity of this point undergo damped oscillations 
and eventually converge to the minimum, which is easily verified by perturbation 
theory. Writing 

(l3{x+) ^ C-^ + S(j){x+) (25) 

and linearizing (|2ip in 6(j) we obtain 



V+d+ - 1 + 2C-2^^^+ S(j) = 



(26) 



The solutions of (l26l) take the form 



Y^ a„e''""=^ 



(27) 



where a„ are arbitrary complex numbers and 



V+s,, = 1 



1 



21nC 



Wn 



41nC 
C2 



(28) 



where Wn are the branches of the Lambert-W function and n runs over all 
integer values. One may verify that the s„ are complex numbers which appear 
in conjugate pairs and also that Re [s„] < for all n. Thus, all solutions which 
approach (f> = C~^ must eventually converge to the vacuum. 

From the preceeding paragraph, we expect that any unbounded solutions 
must be associated with dynamics near the maximum = 0. For example, one 
could imagine starting very close to (/> = with a large velocity that would al- 
low the tachyon to roll over the maximum and down the unbounded side of the 
potential. We have claimed that the condition (1^^ is sufficient to forbid such 
behaviour (and verified this claim also numerically). Let us now seek to under- 
stand better the origin of this constraint. Consider an initial function (po{x^) 
such that 00 (0) is very close to zero. Intuitively, such a function represents the 
most dangerous possibility. By continuity, because (/>o(0) is near the maximum, 
then so is 0i(O). By writing the equation of motion (|2ip in delayed form we can 



^Interestingly, the constraint II23I I is sufficient but not necessary to have sensible evolution 
at late times. 
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establish the following result: 

F+ 9+0x1^+^0 = MO)-C'^li-T) (29) 

> Mo)-M-T) (30) 



where the inequality on the second line follows from (j23l) . If we further assume 
(Pi)) then we have 

5+0iL+=o>O (31) 

so that the solution on the first know must start near the maximum ip — and 
roll towards the true minimum — C~^ . We can describe the early stages of 
this motion using perturbation theory. Writing 

= + d> (32) 

and linearizing (pT|) we find the growing solution 

5(j) = aoe"^/^^ (33) 

The preceeding discussion suggests that our results are rather general. As 
long as (I23I24P are satisfied one expects that all solutions will roll away from 
the maximum, experiencing exponential growth, until they reach the vicinity 
of (/) = C~^ , at which point damped oscillations take over and the solution 
converges to the minimum. Indeed, precisely this behaviour has been observed 
in our numerical analysis. This discussion is also consistent with 28]. 

Finally, let us comment on a "priviledged" class of initial functions. If we 
wish to minimize derivative singularities (which are presumably unphysical), we 
should choose 4)o{x'^) to be an approximate solution of the equation (|2T|) . For 
solutions rolling away from the maximum, we have a one-parameter family of 
such functions: 4>a{x'^) — ee^ 1^ where e ^ 1 quantifies the initial displace- 
ment. With this choice, the condition ([M)) becomes unnecessary, because the 
dynamical equation (|2ip forces solutions near the maximum to grow monoton- 
ically. The possibility to roll down the unbounded side of the potential while 
satisfying (j23|) arises only for initial functions that are very far from being so- 
lutions of (I^TI) . Hence, this possibility will typically also be associated with 
unphysical derivative singularities. 

3.3. Stability and Predictivity 

The dynamics illustrated in Fig. |3] and Fig. IH] are quite remarkable. In both 
cases the Hamiltonian is unbounded from below. Nevertheless, nonlinear so- 
lutions cannot access the unstable directions in phase space starting from an 
"allowed" initial function (j)o(x~^). The constraint which must be imposed on 
4>o{x^) in order to obtain sensible late-time evolution is extremely mild: it 
amounts to restricting the tachyon to the physically meaningful region of the 
potential over an interval comparable to the microscopic scale of nonlocality|j 



''In particularly, we do not require that the field sits at an unstable maximum for an infinite 
amount of time. That would be incompatible with quantum mechanics and, presumably, also 
with macroscopic causality/predictivity. 
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These mild physical restrictions on (j)o{x~^) have the effect of constraining the 
theory to a stable subset of its solution space. That this might be possible has 
been suggested previously, for example by Woodard and Eliezer in [4| . The ap- 
proach employed here has the advantage of being completely nonperturbative 
and also physically well-motivated. Moreover, our solutions depend on infinitely 
many free coefficients and, in this sense, retain information about the full nonlo- 
cal structure of the model. (Hence, our approach is different from perturbative 
localization.) 

Suppose that we had adopted a different view of the initial value problem: 
fixing infinitely many derivatives d\ tp at i = 0, instead of the initial function 
formulation ([6]). In this picture we would have uncovered an "island of stability" 
in initial condition space which leads to stable late-time dynamics. Constraining 
the theory to the island of stability requires a seemingly delicate and contrived 
tuning amongst the infinite tower of initial conditions. However, this tuning 
seems contrived only because one is using a highly unnatural formulation of the 
initial value problem. In light of the initial function formulation ([6|) we see that, 
in fact, the constraints imposed by late-time stability are very mild. 

What becomes of predictivity in our delayed models (O and (pi]) ? As dis- 
cussed in section [21 there is a breakdown of predictivity on microscopic scales. 
Ax"*" ^ T ~ F+. In essence, this occurs because our ignorance about the 
past history of the system makes the initial function (J)q (a:"*" ) inherently ambigu- 
ous. However, as we have seen, the late-time dynamics are largely insensitive 
to (j)o{x'^), provided it belongs to the class of allowed initial functions. Hence, 
we are still able to recover predictivity on macroscopic scales, Aa;"*" 3> T. 

We view the breakdown of microscopic predictivity as a necessary feature of 
working with a fundamentally nonlocal theory. Indeed, it seems quite natural: 
we should not expect to be able to unambiguously resolve dynamics on times 
comparable to the underlying scale of nonlocality. 

Obviously, our conclusions about the nonperturbative stability of equations 
([9|l and (|2ip are purely classical. One might worry that instabilities re-appear 
at the loop level (see [37[ for a related discussion) or that the tachyon field 
(jj can tunnel to the unbounded region of the potential, violating the physical 
restrictions ([TT|) and (|23]) . We note, however, that in string theory the situation 
may be rather subtle. It is not clear if the physical degree of freedom associated 
with the tachyon field even exists on the unbounded region of the potential, or 
at the local mininium|j 

To summarize: the string theory models (0) and (PT|) are nonperturbatively 
stable and make unambiguous predictions over time scales large compared to the 



^At the local minimum the D-brane has disappeared and there can be no open string degrees 
of freedom. It may be that (f> somehow encodes the closed string gas which is the outcome 
of the brane decay |41|| or it may simply be that the </> disappears from the theory. The 
latter interpretation is similar to what happens to the inflaton field during brane/anti-brane 
inflation [ill . At the end of inflation, the mobile brane annihilates with its anti-brane and 
the inflaton field (corresponding to the inter-brane separation) no longer exists as a physical 
state in the theory. 
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delay, provided (j)o (t) belongs to a class of allowed initial function. Obviously, our 
optimistic conclusions have relied heavily on the delayed form of the prototype 
model ([2]) and the formulation © of the initial value problem. Oiily a very 
privileged class of nonlocal theories (those having compact support [4|]) admit 
such a formulation. One might wonder how much of our analysis can be carried 
over to a more realistic QFT framework. We will investigate this question in 
section ID 

4. The Initial Function Formulation Applied to Generic Nonlocal 
Field Theory Models 

Our goal now is to determine the extent to which the optimistic conclusions 
of section [3] carry over to a general nonlocal QFT framework. To this end, we 
first introduce the very general class of nonlocal field theories. Next, we will 
show how the simple-minded inclusion of dissipation into the dynamics leads to 
delay- type nonlocality and an initial value problem of the form ^ . 

4.I. Prototype Nonlocal Field Theory 

Consider the prototype nonlocal field theory 

£=i0F(n)0-^„t (</>), ^„t(0)-|0' + ^</'' + --- (34) 

This model subsumes known examples from SFT [l[, p-adic string theory [6| 
and also cosmology [16|, llSl, [ITI, l2l|, |22|, |23| . The full scalar potential is 






2 



9 

m 






- 2!-^ ■ 3!-^ ■ Ar ■"■•• ^^^^ 

-F(0). 

The nonlocal structure of (p4|) is encoded in the pseudo-differential operator 
F{0). In most cases of interest this can be represented as a convergent series 
expansion 

oo 

F(n) = ^/„n» (36) 

n=0 

where the power D" is to be understood as n-fold composition of the d'Alembertian 
operator with itself. By assumption the generatrix F{z) is therefore an analytic 
function of the complex variable z. The Weierstrass factorization theorem en- 
sures that we can write 

M 



Fiz)^r{z)l[iz~mlY- (37) 



a=l 
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with r(z) having neither zeros nor poles at finite z. The propagator G(j?^ ~ 
l/i^(— p^) therefore has precisely M poles at the locations p^ = —m'^, with the 
a-th pole being of integer order ra- 

The field equation derived from (|M)) takes the form 



F{D)c^^VU^) 



(38) 



The initial value problem associated with (1551) has been studied in [24, |31| | 
using the formal generatrix calculus, and also in [J, [181, |32|, |33|, |3J, [32, 135[ using 
alternative approaches. 



^.^. Multi-Particle Form of the Action 

Let us now re- write the prototype Lagrangian p4p in a form which is equiva- 
lent but which makes certain key features more transparent. (See also 25|, |12|.) 

We assume a kinetic function of the form (157)) and perform the field redefi- 
nitioiH 

^ = r(n)i/20 (39) 

so that the Lagrangian (IM|) becomes 



M 



r-^-^llia^mlr^-v^nt r-i/2(n)^ 



o=l 



Next, we introduce M independent "constituent" fields defined as 



(40) 



(41) 



The "collective" field variable (p of the original theory is a superposition of these 
constituents: 



M 



= r-'/^a)J2va^'' 



(42) 



where the coefficients rja define a partial fraction decomposition of the quantity 
H{z) = F{z)/T{z) as 



Va 



M 

y 

fr{ (z - ml) 



2\ra 



■ M 
.a=l 



-,2\rc 



Now, the Lagrangian p4p can be written in the decoupled form 



M 



M 



r-V2(a)^ 



Vb^ 



6=1 



(43) 



(44) 



"This is an allowed manipulation because r(2) contains neither zeros nor poles. 
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Equation ([^^ is the main result of this sub-section, it shows how the Lagrangian 
([34]) for the "coUective" field cf) can be decomposed into a theory of M "con- 
stituent" fields (^" (one for each pole of the propagator) . Notice that we have not 
made any particular assumptions about the geometry of the background space- 
time, nor have we employed perturbation theory. In Appendix B we show how 
to use the multi-particle form of the action (pij) to construct the perturbative 
Hamiltonian and look for constituent ghost excitations. 

Using (PTj) . P^ and (|33|) is is possible to re- write the dynamical equation 
in the form 



(n~m^r»^, = r-V2(n)v^; 



-1/2 



(°)E 



rib'P 



(45) 



(Notice that here V(^^ is the derivative of l^nt with respect to the collective field 
(f), and not with respect to the constituent (p".) 

4-. 3. Dissipation and Delays 

Our examples in section [3] both inherit their delayed structure from some 
source of dissipation in the dynamics. We can model the effects of dissipation 
on our prototype nonlocal QFT p4p is a very simple way, by adding a friction 
term to the d'Alembertian. Thus, we take 



D = -a^ - pdt 



(46) 



when acting on a homogeneous function. There are many ways to motivate the 
appearance such a friction term. Obviously, we could obtain ()46|) by working 
in an expanding FRW universe. In this case we would have /3 — 3H where 
H = d/a is the Hubble scale. More generally, one might imagine motivating 
the appearance of such a term by supposing that the motion of leads to 
particle production, in which case (3 would be related to the imaginary part of 
the effective action. (See, for example, 48].) 

To see how the simple-minded modelling of dissipation encapsulated in equa- 
tion ()46|) can lead to delay-type nonlocality, let us assume a kinetic function of 
the form ([57]) with 

r(n) = e-"° (47) 

This type of nonlocality is ubiquitous in string theory; factors of the form e~"° 
typically arise at the vertices of Feynman diagrams. Using (|^51) and (|T7)) we can 
re- write the constituent equation of motion (|45p in the form 



{-d^-f3dt-rnlr^^'^{t)^e—^^K 



e-f^?^^,<^&(i-a/3) 



(48) 



where we assume that /? is a constant. We now introduce yet another set of 
basis fieldflj 

^''{t) = e-^^*V'(i) (49) 



'Once again, this manipulation is allowed because e "" '^ has neither zeroes nor poles. 
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Finally, the dynamical equation (P5|) for the constituents takes the form 



i~d^~f5dt-mly-rit) = e 



^ Hit 



J2v'>i'''(t-T) 



where the delay is 



a/3 



(50) 



(51) 



This we suppose is positive. Equation ([50)) is the main result of this sub-section. 
Equation ([50)) is almost exactly of the form ([5]) . The only complication is the 
nonlocal operator e~""« on the right hand side. We can seek solutions by the 
method of steps. (It might be that this approach does not access the full solution 
space.) Each constituent tp"'{t) is split into solution segments ?A'^(i) = ^f (i) on 
the interval t € [(« — 1)T, iT] with integer i. (Here the upper index on the field 
denotes the constituent label while the lower index labels the segment.) The 
j-th solution segment can be obtained recursively as 



{-d^-pdt-mlr^p-{t) = e-^^W^ 



Y,Vb^liit-T) 



L b 



(52) 



On the first segment we have an arbitrary initial function for each constituent 
field 

Vj"(i)=Vo(0 on te[-T,0] (53) 



Since equation (1521) is 2ra-th in time derivatives, it follows that we can make 
ip°'{t) and its first 2ra — 1 derivatives continuous at the knots t — 0,T, 2T, ■ ■ ■ 

It is worth pausing to admire the simplicity of equation ()52p . We can con- 
struct solutions of a nonlinear differential equation of infinite order simply by 
solving a set of M inhomogeneous linear equations, each of which is finite order 
in derivatives. The only complication is the appearance of the nonlocal operator 
e~ * in the source term. 

In this section we have shown that the naive inclusion of dissipation ()46p 
into the general prototype nonlocal theory ([34)) leads to delay-type nonlocality, 
as manifested in the dynamical equation (|50p . Even though the nonlocality is 
not purely of the delay type, we are still able to employ the method of steps. 
In the next section we will show how this approach can be used to construct 
efficient and numerically stable time-dependent solutions. 

Equation ()53p provides a novel new way of thinking about the initial value 
problem for the nonlocal theory ([34)) . Each pole of the propagator contributes 
a "constituent" field (p" with its own independent initial function. 

Of course, the fact that a delayed formulation ([5^ with initial value problem 
([55)) is possible does not imply that the theory ([M)) is stable or usefully predictive 
in any sense. These key questions will depend sensitively on the underlying field 
theory and the class of allowed initial functions. It seems necessary to proceed 
on a case-by-case basis. In the next section, we will consider a specific example. 
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5. A Case Study: Cosmological D-Brane Decay 

5.1. Delayed Formulation and the Method of Steps 

In the last section we showed how the simple-minded inclusion of dissipation 
leads to a delayed formulation of the dynamics in a very general nonlocal QFT 
context. Here we consider a specific example: D-brane decay in a background 
de Sitter space-time. This toy model may be relevant for nonlocal models of the 



early universe [2l|, [2^ |23| and also quintessence [16| . We employ the same level 
zero truncation of SFT as was studied in subsection 13.21 however, this time we 
take a constant dilaton profile and minimally couple (jj to Einstein gravity: 



9s 



3 



i0C-2°(n + i)0-^^ 



(54) 



(see [5] for more details). This example is consistent with our prototype model 
(p4|) and also the more restrictive class of theories considered in subsection 14.31 
The equation of motion 

C-^° (D -t- 1) = C^./)^ (55) 

admits infinitely many initial conditions and the Hamiltonian associated with 
((54|) is unbounded below [24| (see Appendix C for a review). Our goal is to 
derive a delayed formulation of this equation and argue that the theory can be 
constrained to a stable subset of its solution space at the expense of certain mild 
physical restrictions on the class of allowed initial functions. 

In a background de Sitter space-time the equation of motion (I55|) can be 
written as 

{dl+iHdt~l)(t>{t)^-C^'^^^'^{t~T), T = 6HlnC (56) 



Following subsection 14.31 we can employ the method of steps. The i-th solution 
segment 4'i{t) is constructed recursively as 

(92 + 3Hdt - 1) Mt) = -C3"29?02_^ (^ _ J.) (57) 

On the first segment (i — 0) we have the initial function: (j){t) = (t>o(t) on 
t G [—T, 0]. Since the differential operator on the left hand side of ([57)1 is second 
order, we can demand that the field and it's first derivative are both continuous 
at the points t = iT where the solution segments are joined. 

5.2. Allowed Initial Functions 

Which restrictions on 0o(t) are necessary to constraint the theory (|56p to 
a stable subset of its solution space? Following the analysis of subsection 13. 2( 
we might hope that the simple physical requirement < 0o(O < ^ "^ ^ill 
be sufficient. However, this is not the case, we will also need to bound the 
derivatives of (l)o{t). 

We can see which class of (poit) will lead to late-time instability in a very 
simple way, by inspection of the delayed formulation (|57p . For a given initial 
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function (fioit) we can expand the quantity —C^(f)Q{t — T) on the interval — T < 
i <0 as 

- C^(blit -T)=J2 ane^^"'*/^ (58) 

n 

(For the sake of this argument the reader may interpret the summation over 
n as either discrete or continuous.) Now consider the i = 1 solution segment, 
(jji (t) . This obeys an inhomogeneous equation 

{df + 3Hdt - 1) Mt) - E G'"'"'^^'a„e2-"*/^ (59) 

n 

Then the dominant contribution to the source term obviously comes from the 
mode n which maximizes the quantity C^'^ " '"^ a„. Unless the coefficients 
in (j58p tend to zero sufficiently fast, then it may happen that the dominant 
contribution to the source comes from some high frequency mode with n 3> 1. 
If that happens then the particular solution of ([57]) will go like (/)i (t) ^ e*"* with 
uj ^ H . Carrying through the same procedure, the next segment oscillates even 
faster, as (t>2{t) ^ e^*"^*. And so on. At each successive step the energy cascades 
into higher-and-higher frequencies so that the solution will oscillate ever more 
rapidly. This is precisely the usual Ostrogradski instability. It is clear that 
in order to evade this pathological behaviour we must ensure that the Fourier 
coefficients a„ tend to zero sufficiently quickly at high frequency. 

Establishing a necessary condition on (/'o(i) for late-time nonperturbative 
stability is very difficult. Therefore we will instead seek a sufficient condi- 
tion. An obvious candidate is the constraint < 4>o{t) < C~^ and also that 
dl 4>o{t) / (t>o{t) < 7J". In the next section we will test this conjecture using fully 
nonlinear numerical simulations. Although this constraint is slightly stronger 
than those which were considered in section[3l it is still far from onerous. Indeed, 
we might imagine that this type of consistency condition could be established 
dynamically within the context of a more realistic model of dissipation. 

5.3. Numerical Methods and Nonlinear Solutions 

We now implement numerically the method of steps (j57p . Obviously, to 
proceed we need to a faithful numerical implementation of the infinite order 
operator C~^^* . We could proceed by decomposing (j)f_i{t — T) into a Fourier 
series on the segment t € [{i— l)T,iT]. Then each member of the sum is an 
eigenfunction of d^ so that C~^^' acts in a simple way, term-by-term. However, 
this idea turns out to be a very bad approach. The problem is that the Gibbs 
phenomenon introduces errors at the points t — iT where the solution segments 
are joined. These errors propagate and become very large after just a few 
steps. A much better approach is to represent the quantity 4>1_-^{t — T) using an 
expansion in Chebyshev polynomials. When this expansion is truncated to order 
2N (with N integer) then the pseudo-differential operator C^^^* is equivalent 

to the partial sum X]ri=o — t — (^t " ^^'^ i^ straightforward to implement 
numerically. 
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When the dissipation is sufficient to allow the field to settle down to the 
minimum at late times, then our iterative solution of (|56p is very efficient and 
numerically stable, a significant improvement over higher order perturbation 
theory or the diffusion equation formulation (of course, our method is also con- 
siderably less general than those alternative approaches). These results are 
robust when the number of steps or the order of the Chebyshev polynomials is 
increased. On the other hand, when the dissipation is small then the method 
breaks down within just a few steps. However, since we regard wildly unstable 
solutions as unphysical, this is not a serious drawback. 

We have investigated numerically the sufficient conditions on (f>o{t) for late- 
time stability. The results confirm our previous intuition. We will have </> — >■ C~^ 
as t — > oo as long as < 0o(i) < C~^ on t e [— T, 0] and the derivatives of (j)o{t) 
are not too large as compare to the coefficient of dissipation!^ These dynamics 
are illustrated in Fig. [T] 



Time -Like Cosmologlcal Tachyon with H =2 m^. 




Figure 7: The time-evolution of the time-like SFT tachyon in the presence of cosmo- 
logical (de Sitter) friction for several representative choice of initial function, c/>o(t). 
We set H = 2ms for illustration. Provided the constraint < 00 (i) < C~^ is satisfied, 
and also that the derivatives of <j)o{t) are not too large, then the field will eventually 
settle down to the true minimum at late times. The horizontal yellow lines illustrate 
the critical points (j) — and (j) — C~^. The explicit choices of 4>o{x'^) are provided in 
Appendix A. 



Our results suggest that nearly any source of dissipation will allow us to 
salvage sensible dynamics from the theory (j54p . even without making restric- 
tive assumptions about the symmetries of the solution. We find this prospect 
encouraging since the rolling tachyon is supposed to be decaying into closed 
string excitations [40|| which were not taken into account in the effective de- 



^^ Quantitatively we can demand that the Laplace transform of (/>o(t) is analytic inside a 
circular contour of radius H. We stress that this is a sufficient, but not necessary, constraint. 
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scription ([M|) . Thus, our results may provide a plausible mechanism by which 
the consistent dynamical inclusion of closed strings could allow SFT to evade the 
Ostrogradski instability. At the more practical level, this general prescription 
helps to establish the viability of nonlocal cosmological models. 

5.4- Wild Oscillations and Allowed Initial Functions 

It has been argued previously that the dynamics of tachyon condensation are 
stable for th e sp ecial case where (j>{t) — >■ as i — >■ — oo. For example, this was 
discussed in [29| for the light-like dilaton models discussed in section |3] (although 
in that paper the analysis went beyond our simple level-0 truncation). See also 



32|. It is not clear if demanding that the field sit at an unstable maximum 
can be consistent with quantum mechanics (or, for that matter, macroscopic 
predictivity) . It should therefore be emphasized that our class of "allowed" 
initial functions is much broader than this very special subset of solution space. 
Wc only require certain qualitative restrictions on the dynamics of (f>{t) over a 
microscopic interval At ^ TnJ^. Moreover, our solutions depend on infinitely 
many free parameters and, in this sense, encode information about the full 
nonlocal structure of the equations. Nevertheless, it is worth comparing our 
analysis with the special case where (j){t) — >■ as t ^- — cx). Because = is a 
surface of degeneracy, there is only a single growing mode which can be accessed 
by naive perturbation theory about the unstable maximum: (f>Q{t) — eC~^e^+*, 
with e <C 1 and r+ = — ^ + ^^9H^ + 4. Using the method of steps, we have 
confirmed numerically that this special class of initial function always leads to 
sensible late-time evolution, at least when i? > 1. Notice that this special 
function is not the unique perturbative solution of equation (j56p which satisfies 
4>{t) — ?► as t — ^ — oo. Using the technique described in the appendix of [28| 
we could easily constuct a family of such solutions which depend on infinitely 
many free parameters. 

In passing, it is worth noting that the wild oscillations of rolling tachyon 
solutions, which are a symptom of the Ostrogradsi instability, are not always 
viewed as pathological. This comes down to a delicate issue of interpretation. 
It is nontrivial to relate 4> to physical observables in SFT. Perhaps the ghostly 
constituent excitations are pure gauge. In ^26| it was argued that the unstable 
dynamics of can be interpreted as describing the closed string physics asso- 
ciated with Sen's tachyon matter J49|. It is not clear to us how (or if) this 
claim addresses the issue of predictivity discussed in the introduction. (See also 

50, l5l|, l52] for related discussions.) In this work our aim is to address stabil- 



ity and predictivity in a general QFT context, without relying on intrinsically 
stringy physics. Hence we regard any wildly unstable solutions as unphysical. 

6. Conclusions 

In this work we proposed a new formulation of the initial value problem for 
nonlocal theories with infinitely many derivatives. This new approach is less 
general than the pseudo-differential operator technique advanced in [2J, Isij . 
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however, it has considerable conceptual appeal. In particular, the underlying 
nonlocality of the model is manifest; the Cauchy surface has been "smeared 
out" over a finite width. 

We have found that this mathematical re-formulation leads to a novel new 
perspective on the key problems of stability and predictivity which plague non- 
local models. By imposing mild physical restrictions on 4>a{t)^ it is possible to 
constraint the theory to a stable subset of its solution space. Moreover, the con- 
strained theory makes robust predictions for dynamics on macroscopic scales. 
Although our claims are limited only to a small number of specific examples, 
we expect that similar results might be possible more generally. It is our view 
that these restrictions on the initial function should be viewed as a part of the 
definition of the theory. 

Philosophically, our prescription is very similar to approaches which have 
been advocated previously for higher derivative theories of finite order. Both 
the Cutkowski et al. contour deformation 43|, |4J, |45| and also the Hertog and 



Hawking [46| approach are, classically, equivalent to adopting stability as a 
future boundary condition. Both these approaches lead to a violation of mi- 
croscopic causality, however, it is thought that macroscopic causality survives 
46l . |47| . In the initial function formulation we have an elegant and physically 



well-motivated way to impose late-time stability for models with infinitely many 
derivatives. Our approach is nonperturbative and retains information about the 
full nonlocal structure of the theory. Notice that microscopic causality is man- 
ifestly violated in this formulation of the initial value problem. We view this, 
along with the breakdown of microscopic predictivity, as a necessary feature 
of working with a fundamentally nonlocal theory. In fact, these features seem 
rather natural: we should not expect to be able to unambiguously resolve dy- 
namics on times comparable to the underlying scale of nonlocality. 

It is not our goal to claim that nonlocal models are generically sensible or 
phenomenologically viable. Rather, we believe that fundamental nonlocality 
is ubiquitous in well-motivated microscopic models and, therefore, it behooves 
us to try to understand such constructions from a general perspective. By 
studying the connection between dissipation and delays, we have arrived at a 
new formulation of the initial value problem and a novel new perspective on the 
key issues of stability and predictivity. These general lessons may be applied in a 
variety of contexts. For example, the rolling open string tachyon is, physically, 
supposed to be decaying into closed string excitations [40| which are usually 
not taken into account in SFT. Plausibly, the dissipation induced by this type 
of particle production has a dynamical effect which is similar to the dilaton 
gradient in section [3] or the cosmological friction in section [5j Thus, we have 
elucidated a mechanism by which the consistent inclusion of closed strings might 
enable SFT to escape the Ostrogradski instability. A strength of our approach 
is that we do not rely on any intrinsically stringy features of the underlying 
model. Thus, our results may readily be applied to more general nonlocal QFT 
models (in particular, to nonlocal cosmological models). 

There are a number of directions for future work on this subject. Obvi- 
ously, our results for nonlinear stability have been rather specific and niodel- 
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dependent. It would be very interesting to understand the connection between 
dissipation and delays in more general models, for example where the coefficient 
of dissipation is time-dependent (as would arise in a self-consistent cosmology). 

Our analysis has been entirely classical, which is an obvious limitation. In- 
stabilities might re-appear at the loop level, or it might be that the field can 
tunnel quantum mechanically to some unstable configuration even when the 
unstable directions in phase space are not classically accessible from an allowed 
initial state (perhaps similarly to [53|). At the current level of understanding, 
it seems that perturbation theory is best suited to study quantization while the 
method of steps is best suited to study nonperturbative dynamics in the classical 
regime. However, it might be possible to incorporate quantization directly into 
the delayed formulation. One could imagine imposing (j)o{t,y:),(j)o{t,x) 7^ on 
the initial segment — T < ^ < and then using the method of steps to construct 
the q-numbcr field (j){t, x) on i > 0. It would be interesting to understand how 
the imposition of physical constraints on (/)o can be performed in a quantum 
theory. 

One might also imagine trying to compute statistical properties of the t > 
dynamics by performing a functional average of the solution over allowed 0o(O- 
We find this idea fascinating, however, it is not clear how to put a measure on the 
space of initial functions. We have investigated this idea using the exact solution 
(|12p where the averaging can be performed analytically. In this particular case 
it seems that the ambiguity in choosing the averaging procedure is not less sever 
than the ambiguity in choosing a realization of the initial function. 

It is interesting to note that our physical restrictions on 0o(i) are not so 
stringent that they forbid violations of the Null Energy Condition (NEC). In 
[28] this was shown explicitly for the toy p-adic model studied in section [3l 
Violations of the NEC play a key role in cyclic/bouncing cosmologies, such 
as the new ekpyrotic model [5J], but are notoriously difficult to obtain in a 
consistent microscopic framework. It is possible that our formalism could be 
used to construct sensible, NEC-violating nonlocal field theories. We leave this, 
and other interesting possibilities, to future investigations. 
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APPENDIX A: Initial Functions for Numerical Analysis 

In this appendix we summarize the explicit initial functions used to generate 
the numerical results Fig. |6] and Fig. [71 In the left panel of Fig. [6] the solid red, 
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dotted black and dashed green curves correspond, respectively, to the functions: 



,{x+) = 10- 






0.25 e^"+^ 
0.25e-3(-++^) 



In the right panel of Fig. [5] the solid red, dotted black and dashed green curves 
correspond, respectively, to the functions: 



M^^) = 0.1 

0o(x+) = -1 
,/,o(x+) = 10^2 sin 



1-lOe 

"27rx+ 



25{x++Tf 



T 



In Fig. [7] the solid red, dotted black and dashed green curves correspond, re- 
spectively, to the functions: 



0.25 6+'^+* 

1 + 0.85t2 cos (i/2) e*/2 



-3gr+t 



Mi) - 10"' 
00 W = 10 
where r+ = — ^ + i\/9-ff^ + 4, as defined in subsection 15.41 

APPENDIX B: Conditions for a Bounded Hamiltonian 

In general the interactions in (HH) will be nonlocal and the dynamics are 
very complicated. We can gain some insight into the dynamics of the theory by 
adopting a perturbative approach. We write the each constituent field as 

¥>" = Sif" (B-1) 

and expand the Lagrangian ()44p order-by-order in perturbation theory 

C = (5£(2) + 5£(3) + . . . (B-2) 

where the quadratic part 50-'^' gives the linearized equations of motion and 
the higher order contributions encode the effects of interactions. The quadratic 
Lagrangian decomposes into constituents as 

j\/ 
5C^'^ = Y: 5Cf\ SCi'^ = ^S^^in - mly-^S^^ (B-3) 

a=l 

Obviously the quadratic part of the Hamiltonian must also admit a similar 
decomposition 



M 



dn'' 



E^^i' 



(B-4) 



o=l 
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where each SWa is constructed from SCa ' . 

Under what conditions can the perturbative Hamihonian ()B-4|) Hamihonian 

(2) 

be bounded below? If r^ > 1 for any a then 5Ca must depend non-degenerately 
on time derivatives higher than than first order. By the Ostrogradski theorem 

(2) 

it fohows that SHa must be unbounded below. Hence, only theories with 



Ta = 1 for all a yield a bounded Hamihonian. In this case the coefficient ria can 
be computed easily from the residues of the propagator at the pole m^: 






(The contour Ca only encloses the pole at z = m.^.) If any rja are negative, then 
((311) must contain ghosts. If Af > 1 then, because H{z) is analytic, the residues 
must flip sign and some T-La will be unbounded below. Therefore the only pos- 
sibility to construct a perturbatively ghost-free theory is when the propagator 
G(p^) ^ l/F{—p^) contains at most a single pole. 

We have implicitly assumed that rn^ is real valued for all a. The case of 
complex mass-squared is slightly more subtle since now the rja coefficients are 
complex, as are the fields Sf"' . Hermiticity of the action requires that such 
complex poles arise in conjugate pairs, so that if rn? is a zero of F{z) than so 
must be {m?)* . Let us suppose that m^ and rn^ are complex conjugate pairs. 
Then, reality of the solution implies that Sip^ — {Si^^)* and the relevant part of 
the Lagrangian is 

'^5^\U - ml)5^^ + c.c. = -]^KABd^,5^k''d^5i>'' - ^^MabH^H"" (B-6) 

where c.c. denotes the complex conjugate of the preceding term and in the 
second equality we have defined S^pi = -h= [Sijj^ + iSijP'^ so that the indices 
A, B run over real and imaginary parts. The kinetic term can be diagonalized 
by a rotation 5^^ — >■ RgS^p^ . It is easy to verify that one eigenvalues of Kab 
is always negative and the Hamiltonian is unbounded below. 

To summarize the results of this appendix; we have shown that a necessary 
(but not sufficient) condition for ([M)) to yield a perturbatively stable Hamilto- 
nian is that the propagator G{p^) ~ 1/ F{—p^) has at most a single pole. In 
this case the theory describes only a single degree of freedom and the equation 
of motion (p8| requires only two initial dataj^ Absent the imposition of some 



constraints (for example on the initial function 0o(O o^ following |43l.l44ll45l.l46j) 
ghost-like constituent excitations will lead to catastrophic classical instabilities 
and vacuum decay in the quantum theory. See [42| for stringent observational 
constraints on the existence of physical ghost fields. 

There is one scenario where the analysis of this appendix can give very 
misleading results: the case where = is a surface of degeneracy for the theory 



^^ Stable theories describing multiple degrees of freedom can be constructed when the gen- 
eratrix is not analytic, for example if F{z) has isolated poles or branch cuts. 
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([M]). In this case the propagator obtained in a naive perturbation theory about 
(/> = will have less poles than the total number of physical degrees of freedom 
in the theory. This situation actually occurs for the level truncation of SFT, see 
appendix B for a review. 

APPENDIX C: Counting Initial Data for the SFT Tachyon 

It is sometime incorrectly argued that (j55p admits only two initial conditions. 
This specious conclusion is reached by perturbing about the constant solution 
(j) = 0, which corresponds to the unstable maximum of the potential. Writing 
(f) = + S(j> we have 

C-2°(n + l)(50 = O (C-1) 

to linear order in Scj). The generatrix is F{z) — C^'^^^{z + 1) which has only a 
single zero corresponding to a tachyonic excitation with mass-squared Af ^ = — 1 
(in string units). At higher order in perturbation theory the counting of initial 
data is the same. However, this conclusion is very misleading. Notice that 
perturbing about the other critical point gives a different answer. Writing 
(f) = C~^ + 54> E^iid linearizing \n 54> we obtain 

[C-^° (D + 1) - 2] ^0 = (C-2) 

The generatrix is now F{z) — C^'^'^{z + 1) — 2 which has infinitely many zeroes 
F{Ml) = given by 



Mi 



' .w. ( ''-'' 



21nC V C2 



(C-3) 



where n = 0, ±1,±2, • • • and Wn{x) denotes the branches of the Lambert-W 
function. The spectrum (jC-3|) corresponds to an infinite tower of constituent 
physical states with complex mass-squared. Correspondingly, equation (jC-2p 
admits infinitely many initial conditions. In the Appendix of [28| we showed 
how to construct solutions of ()55|) near = which depend on infinitely many 
free coefficients. Using the results of Appendix A, it follows that the Hamiltonian 
associated with ([51]) must be unbounded below. 

Perturbing about different critical points of the potential we found different 
numbers of initial conditions. This is because = is a surface of degeneracy 
for the equation of motion (j55p 4] . This phenomenon has nothing to do with 
the nonlocal structure of the theorvr^l 



^^A simple example of a local theory with the same property was provided in 4]. Consider 
the Lagrangian £ = — ^^-'lp■ip'^ . The Euler-Lagrange equation is fourth order in time deriva- 
tives and hence generic solutions admit four initial data. However, on the surface ?/) = the 
character of the equation of motion changes: instead of determining d'^ip in terms of lower 
order derivatives it instead becomes a constraint on the lower derivatives. Thus, solutions 
which start at the point i/i = contain only three free coefficients. 
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The fact that (|51)) allows for negative kinetic energy cannot be doubted. 
Time-dependent solutions in flat space have been shown to move higher on 
the potential than their starting point, without difficulty 7]. Since energy is 
conserved in this rolling process, it follows that the kinetic energy must be 
negative. Moreover, generic solution of ([55]) undergo erratic oscillations with 
divergent frequency and amplitude, the hallmark behaviour of an Ostrogradski- 
sick theory. The erratic time dependence of rolling tachyon solutions is properly 
understood as the symptom of a much deeper problem with SFT. 
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